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Abstract 

We determine the 2-loop effective gauge coupling of QCD at high temperatures, defined as a 
matching coefficient appearing in the dimensionally reduced effective held theory. The result 
allows to improve on one of the classic non-perturbative probes for the convergence of the 
weak-coupling expansion at high temperatures, the comparison of full and effective theory 
determinations of an observable called the spatial string tension. We find surprisingly good 
agreement almost down to the critical temperature of the deconfinement phase transition. 
We also determine one new contribution of order O(g^T^) to the pressure of hot QCD. 


March 2005 



1. Introduction 


Indirect signs for rapid thermalisation after heavy ion collisions at RHIC energies, derived for 
instance from the fact that hydrodynamic models assuming local thermodynamic equilibrium 
appear to work very well 0, have underlined the need to understand the physics of thermal 
QCD at temperatures above a few hundred MeV. 

Given asymptotic freedom, a natural tool for these studies is the weak-coupling expan¬ 
sion j2j. Alas, it has been known since a long time that the weak-coupling expansion con¬ 
verges very slowly at all realistic temperatures Ham. It also has theoretically a non-trivial 
structure, with odd powers of the gauge coupling .5] and even coefficients that can only be 
determined non-perturbatively 00 . 

On the other hand, the degrees of freedom responsible for the slow convergence can be 
identified umm- they are the “soft” static colour-electric modes, parametrically p gT 
(leading to the odd powers in the gauge coupling), as well as the “ultrasoft” static colour- 
magnetic modes, parametrically p ~ g 2 T (leading to the non-perturbative coefficients in the 
weak-coupling expansion). Here p denotes the characteristic momentum scale, g the gauge 
coupling and T the temperature. The belief has been that perturbation theory restricted 
to parametrically hard scales p ~ 2vrT alone should converge well, while the soft and the 
ultrasoft scales need to be treated either with “improved” analytic schemes, or then non- 
perturbatively. As a starting point for these demanding tasks one may take, however, either 
the dimensionally reduced effective field theory mm or the hard thermal loop effective 
theory m, which have been obtained by integrating out the parametrically hard scales. 

Quantitative evidence for this picture can be obtained by choosing simple observables which 
can be determined reliably both with four-dimensional (4d) lattice simulations and with the 
soft/ultrasoft effective theory. This forces us to restrict to static observables and, for the 
moment, mostly pure gauge theory. Various comparisons of this kind are summarised in 
Refs. fTH m ITT)] . The most precise results are related to static correlation lengths in various 
quantum number channels EH, where good agreement has generally been found down to 
T ~ 2T c , where T c is the critical temperature of the deconfinement phase transition. The 
thermodynamic pressure of QCD is also consistent with this picture EH, even though that 
comparison is not unambiguous yet, due to the fact that the effective theory approach does not 
directly produce the physical number, but requires not-yet-determined ultraviolet matching 
coefficients for its interpretation m 1 

The purpose of this paper is to study another observable for which an unambiguous compar¬ 
ison is possible. The observable is the “spatial string tension”, cr s . 4d lattice determinations 
of a s in pure SU(3) gauge theory exist since a while already [22J but, as has most recently 
been stressed in Ref. 1231 , the comparison with effective theory results shows a clear discrep¬ 
ancy. In order to improve on the resolution on the effective theory side, we compute here 
the gauge coupling of the dimensionally reduced theory up to 2-loop order. Combining with 

1 For the status regarding a few other observables, see Refs. EMEU. 
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other ingredients EHHSI, to be specified below, allows then for a precise comparison. We find 
that once the 2-loop corrections are included, the match to 4d lattice data improves quite 
significantly and supports the picture outlined above. 

The plan of this paper is the following. In Sec. [2] we present the 2-loop computation of 
the effective gauge coupling of the dimensionally reduced theory. In Sec. El we discuss the 
numerical evaluation of this result. In Sec. 0 we use the outcome for estimating the spatial 
string tension, and compare with 4d lattice data. We conclude in Sec. El 


2. Effective gauge coupling 

We consider finite temperature QCD with the gauge group SU(IV C ), and Nf flavours of mass¬ 
less quarks. In dimensional regularisation the bare Euclidean Lagrangian reads, before gauge 
fixing, 

Sqcd = f 0 ^ T f d^^QCD, (2.1) 

£qcd = ^ F ^ F pu + ihiLDvrt > ( 2 - 2 ) 

where (3 = T _1 , d = 3-2e, n,v = 0,..., d, F* v = Dp = d^-igpA^, 

Afj = A^T 11 , T a are Hermitean generators of SU(IV C ) normalised such that TV [T a T b ] = 5 ab / 2, 
7 ^ = 7 n, { 7 ^, 7 ^} = 2 6nu, gs is the bare gauge coupling, and carries Dirac, colour, and 
flavour indices. We use the standard symbols Ca = N c , Cf = (N 2 — 1)/(2IV C ), Tp = Nf/2 for 
the various group theory factors emerging. 

At high enough temperatures, the dynamics of Eq. (EH1) is contained in a simpler, dimen¬ 
sionally reduced effective field theory mm 3 Ei: 

'S'eqcd = J d d x£EQCD, (2-3) 

£eqcd = ^ F ij F ij + Tr [Di, Bo ] 2 + m|TV [Bq] + A^(Tr [Bq ]) 2 + A^Tr [Bq] + ... . (2.4) 

Here i = 1,..., d, Fij = diBj — djB f + gpf abc B b Bj , and Di = di — ig-pBi. The fields have 

the dimension [GeV] 1 / 2-6 , due to a trivial rescaling with T 1 / 2 . Note also that the quartic 
couplings Ap 2 "* are linearly dependent for N c < 3, since then Tr [Rq] = ^ (Tr [Bq]) 2 . 

The theory in Eq. E3> has been truncated to be super-renormalisable; that is, higher order 
operators m (see also Refs. eeiimi and references therein) have been dropped. The relative 
error thus induced has been discussed for generic Green’s functions in Ref. ED], and for the 
particular case of the pressure of hot QCD in Ref. M- In the following we concentrate on 
an observable dynamically determined by the colour-magnetic scale p ~ g 2 T, and it is easy 
to see that in this case the higher order operators do not play any role at the order we are 
working. 
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The effective parameters in Eq. dHH) can be determined by matching, that is, by requiring 
that QCD and EQCD produce the same results, within the domain of validity of the latter 
theory. It is essential that infrared (IR) physics is treated in the same way in both theories at 
the matching stage and, as outlined in Ref. 0, the most convenient implementation of this 
requirement is to perform computations on both sides using “unresummed” propagators. We 
follow this procedure here. 

The matching simplifies further by using the background field gauge (Ref. and refer¬ 
ences therein). As this is essential for what follows, we start by briefly recalling the basic 
advantage of this approach. For a concise yet rigorous overview of the technique, see Ref. m- 
We denote the background gauge potential with and the gauge-invariant combination 
following from symbolically as B 2 + gB 3 + g 2 B 4 . Now, the computation of 

the effective Lagrangian by integrating out the hard scales p ~ 2-kT produces, in general, an 
expression of the type 

-Ceff ~ C 2 B 2 + C 3 gB 3 + C 4 g^B 4 + ... , (2-5) 

where q are coefficients of the form q = 1 + 0(g 2 ). As the next step we are free to define 
a canonically normalised effective field B e g- as B 2 S = C 2 B 2 . Then the effective Lagrangian 
obtains the form 

£eff ~ B 2 S + C 3 C 2 // gB 3 s + C 4 C 2 2 g 2 B 4 s + ... . (2-6) 

We can now read off the effective gauge coupling from the gauge-invariant structure: 

- 3/2 1/2 -1 /o 

9eS = C 3 C 2 9 = c 4 c 2 g . (2.7) 


We observe that two independent computations are needed for the determination of g e g, but 
we can choose whether to go through the 3-point or the 4-point function, in addition to the 
2 -point function (that is, using C 3 or c 4 , in addition to C 2 ). 

The background field gauge economises this setup. Indeed, the effective action is then 
gauge-invariant not only in terms of B e s, but also in terms of the original field B EH- 
Writing Eq. (12.51) as 


-Ceff ~ C 2 


B 2 + C 3 C 2 1 gB 3 + c 4 c 2 1 g 2 B 4 


+ •• 


( 2 . 8 ) 


gauge invariance in terms of B now tells us that C 3 
Eq. (12.71) . we obtain 


_ -1/2 
SeS c 2 g , 


C 2 and c 4 


C 2 - Combining with 


(2.9) 


so that it is enough to carry out one single 2-point computation, in order to obtain g e g. In 
our case, the role of < 7 e ff is played by (cf. Eq. (12.41) 1. 

The class of background held gauges still allows for a general (bare) gauge parameter, £. 
As a cross-check we have carried out all computations with a general £, and verified that it 
cancels at the end. To be definite, we denote (£)here = 1 — (£)standard, so that the gauge held 
propagator reads 





, q^qu 

V) 2 - ' 


( 2 . 10 ) 
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Figure 1: The 1-loop and 2-loop self-energy diagrams in the background field gauge. Wavy 
lines represent gauge fields, dotted lines ghosts, and solid lines fermions. The 2-loop graphs 
have been divided into two-particle-irreducible and two-particle-reducible contributions. 


In order to match the effective gauge coupling, we need to compute the 2-loop gluon self¬ 
energy, II i_ w (p), for the background gauge potential B “. The graphs entering are shown in 
Fig. in The external momentum p is taken purely spatial, p = (0,p), while the heat bath is 
timelike, with Euclidean four-velocity u = (1,0), so that u ■ u = 1, u ■ p = 0. In this case 
has three independent components (Iloi, IRo vanish identically), 

noo(p) = n E ( P 2 ), iiij(p) = (s^ - ^)n T ( P 2 ) + ^n L ( P 2 ), ( 2 . 11 ) 

where i , j = 1, ...,d. In fact loop corrections to the spatially longitudinal part IIl also vanish, 
so that only two non-trivial functions, IIe^t, remain. 

Since we are carrying out a matching computation, any possible IR divergences cancel as we 
subtract the contribution of EQCD. Therefore we may Taylor-expand n /(i , (p) to second order 
in p 2 . This leads to the nice simplification that the results on the EQCD side vanish identically 
in dimensional regularization, due to the absence of any mass scales in the propagators. Thus 
we only need to compute unresummed integrals on the QCD side. 

After the Taylor-expansion, the 2-loop QCD integrals can all be cast in the form 


I(h,i2\ji,32 , j'd ; h,k 2 ,k 3 ) = 


9o r o (q-P) J1 (r- p) 22 (q- r ) 23 


q,r 


[Qo + q 2 ] fcl [ r 0 + r2 D 2 Kvo + r 0 ) 2 + (q + r) 2 ] fc 3 


( 2 . 12 ) 


The indices here are non-negative integers, and the measure is the standard Matsubara sum- 
integral (bosonic or fermionic), with the spatial part /d rf q/(27r) d fd d r/(2ir) d . 
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To reduce integrals of the type in Eq. to a small set of “master integrals”, we 

employ symmetries following from exchanges of integration variables, as well as general partial 
integration identities for the spatial parts of the momentum integrations. The implementation 
of these identities follows the procedure outlined by Laporta m, in analogy with Ref. M- 
We are lead both to very simple 1-loop recursion relations, such as 

/(2*r, 0; 0,0, 0; h ,1,0) = 2kl ~ 2 ~ d I(2i 1 - 2, 0; 0, 0,0; k x - 1,1, 0) , (2.13) 

Zk\ — z 

as well as well-known but less obvious 2-loop ones |2H, like 

7(0,0; 0,0,0;l b ,l b ,l b ) =0, (2.14) 


where the subscripts refer to bosonic four-momenta. 

After this reduction, only six master integrals remain: 

4(n) = ’ I,in) = twr' (2 ' 15) 

where %, qf refer to bosonic and fermionic Matsubara momenta, respectively, and n = 1,2,3. 
For a vanishing quark chemical potential, as we assume to be the case here, the fermionic 
integrals reduce further to the bosonic ones, 


h{n) = (2 2n ~ d - l)/ b (n) , 

leaving only three master integrals. They are known explicitly, 


4 («) 


27T d /2 T 1 + d T(n - d/2) 

(2 TTT) 2n r>) C(2n 


d). 


(2.16) 


(2.17) 


This expression is easily expanded in e and, in the following, we need terms up to 0(e). For 
completeness, the relevant expansions are shown in Appendix A. 

Writing now the Taylor-expanded bare 2-point function IIt of Eq. (EH) as 


IIt(p 2 ) = IIt(0) + p 2 n T (0) + ... 

oo oo 

= E n Tn(0)(9|) n + P 2 E n Tn(0)(5|) n + - , (2-18) 

n =1 n =1 

where gs is the bare gauge coupling, and correspondingly for IIe, our results read 


IIti(0) 

n^(o) 

IIt2(0) 

n T2(0) 


0, 

^-E^/ b ( 2 ) + ^T F 7 f (2), 

b 3 

0, 


(d — 3)(d — 4) 
d(d-2)(d-5)(d-7) 



2ld-7)C\ll(2)~ 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 
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4 C F + (< d 2 - 6d + l)CUjT F J b (2)J f (2) - 
(d 3 - 12d 2 + 39d - 12 )C a - 
- 2(d 3 - 12d 2 + 41 d - U)C f ]t f I?(2)X + 


+ 3^d-7) < i ^ ~ 3W + U4) 4TF/f(1) - ( d “ CU4(3) - 


n E i(o) 

n Ei(0) 

n E 2(o) 


- 8(d - 1 )(d - 6)CyZ> [l b (l) - I f (l) I { (3) 

-(d-l)[4T F I f (l) - (d-l)C A I b (l) 

'd 2 -5d + 28 


6 


+ (d — 3)£ 


CU 4 ( 2 ) + 2(d 0 — r F j f ( 2 ), 


(d - l)(d - 3)^ (1 + 0 [4Tp/ f (l) - (d - 1)CU4(1) CU4(2) + 
+ 4C F r F [4(1) - 4(i)l i f (2) 


( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 


We leave out the lengthy expression for ng 2 (0), as it is not needed in the following. 

The bare results need still to be renormalised. The bare gauge coupling is written as 
g 2 B = g 2 (g)Z gi where g 2 (g) is the renormalised gauge coupling, g is an MS scheme scale pa¬ 
rameter introduced through g 2 = /i 2 e 7 E/47r, and the combination g~ 2e g 2 (g) is dimensionless. 
Denoting 


the factor Z g reads 


Po = 


Pi = 


1 Po -2 


-22C A + 8 T f 


—68 C 2 a + 40 C a Tf + 24 CfTp 


Zq — 1 + — — 2e 9 2 {^) + 


1 


Pi , m: 

4e 4e 2 


(47r) 2 2e r (47t) 4 

and the renormalised gauge coupling satisfies, in the limit e —> 0, 


g-^g\g) + 0{g b ) 


= TpT^9 4 (g) + 7 TTlS ,6 (^) + 0( y9 S ) 


dg 


(47T 


(47r) 


(2.26) 

(2.27) 

(2.28) 

(2.29) 


To proceed, we first cross-check our results for II E against known expressions. After 
the fields Bq of EQCD are normalised to their canonical form (cf. Eq. (12.61) 1. (£>q.Bq) e = 
(BSBZ) 4 d[l + Ilg 1 (0)]/T, we obtain for the matching coefficient m|, 


m E — sb n E i(o) + g% n E2 (o)-nEi(o)n E i(o) + o(g B ). 


(2.30) 
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Inserting Eqs. CSZ5MZZ5) , the ^-dependence duly cancels. Re-expanding also g 2 B in terms 
of the renornralised gauge coupling, and writing then m 


m| = 2*^(0) 


aE4 + Oi E5£ 


+ 


ff 4 (A) 

(47t) 2 


«E6 + /?E2e + C , (s b ! Of’ 


(2.31) 


we recover the known values of aE 4 , «E 5 and «e 6 0) (for original derivations, see Ref. jHJ 
and references therein). We also obtain 


1 


/?E 2 = — Cl 264 In 2 


| 2 641n 2 ^ 




+ 


80 - 1767e + 176 


C'(-i) 


36" "l V 47 tT 

+8 + ll7r z — 887 ^ — 407e — 17671 + 40 


+ CfTf{ — 


{- 81n (S?) _2 “f ln2 + 47E “ 4 


C(-1) 

C'(-i) 

C(-1) 

C'(-i) 

C(-1) 


In 


4ttT J 


+ 


+ 


+ 


+ -OTjl681n 2 (^' 
36 l V47 tT 


+ 


232 - 432 In 2 - 112 te + 112 


C ; (-i) 

C(-i)J 


In 


^e 7 E \ 

4t tT ) 


+ 


+28 + 7t r 2 + 24 In 2 - 64 In 2 2 - 567 ! - 72 7e + I 287 E In 2 - II 271 + 

C'(-l) C'(-l) 

+72 7 ; . / — 128 In 2^ ^ j 


+ 


T%{ -24 In 


C(-i) 

2 / +e 7 E \ 


C(-i) 


+ 


\ 4vrT J 


+ 


8-48 In 2 + I 67 E - 16 


C'(-i) 

C(-i) 


In 


/re 7 E \ 
4vrT / 


+ 


+4 — vr 2 — 8 In 2 + 16 In 2 2 + 87 I — 87 E + 327 e In 2 + I 671 + 
+ 8 ^.f~V - 32In2^ 


C(-i) 


C(-i) 


(2.32) 


Here 71 is a Stieltjes constant, defined through the series £(s) = l/(s — l) + X)Oo7n( — l) n (s — 
l) n /n\. (Note that the Euler gamma-constant is 7 e = 70 -) The result in Eq. (12.321) . first 
obtained in Ref. m by employing the results of Ref. |5j, contributes to the pressure of hot 
QCD at 0(g e T 4 ) [TOj. We rewrite the expression here, since Ref. EH employed an extremely 
compactified notation. 

We then move to consider the transverse spatial part, nx(p 2 ). According to Eq. (12.01) . this 
directly determines the effective gauge coupling: 


5e — T — 9b n(ri (0) + g% (Uti 


- n: 


T2(0) + 0(g B )} . (2.33) 


Re-expanding again in terms of g 2 (g), we parameterise the result (following Ref. jTUJ) as 


9e — Tig 2 (g) + 


2/-\ 


0+) \ 


(47t) 2 


«E7 + P E3£ + 0(e 2 ) 


+ 


g 6 (+) 

(47t) 4 


7ei + 0(e) +0(g 8 )\ . (2.34) 


We recover the known expression H3H2! for q:e 7 , 


«E7 = —A) In 


Bfl) + I C *- 1 *T f In 2 
47 tT 7 3 3 


(2.35) 
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and obtain the new contributions 


/?E3 = 


1 

12 

1 

3 


C A 

T f 


88 \ n 2 (^^) + 8In( ) + llvr z - 887I - 17671 


ge 1 ^ \ 
4 ttT ) 


81n z f + 32 In 2 In f + vr 2 + 16 In 2 2 - 87 ^ - I 671 


7ei = 


M&b 


4vrT J 

21 ge^ \ 

\4irT J ' \AttT J 

fge^\ 1 
3 


(2.36) 


^{^f)-^ C A + ^ T Fln2 


16. 

y 


1 


- —i C A 

18 ' A 


2 r -341 + 20C(3) + 4 C a Tf [43 + 24 In 2 + 5C(3) 
+6 C f T f [23 + 80 In 2 - 14C(3) 


+ 


(2.37) 


The first one, /?E 3 , constitutes again an 0(g 6 T 4 ) contribution to the pressure of hot QCD [TDl . 
while the latter one is the desired finite 2-loop correction to the effective gauge coupling. 


3. Numerical evaluation 


We wish to compare numerically the 1-loop and 2-loop expressions for g Pj , in the limit e —> 0. 
When carrying out such a comparison, it is important to specify the definitions of the 
parameters. Following standard procedures, we solve Eq. (1+291) exactly at 2-loop level, and 
define 


A-ms = Jim g bog 2 {n) 

/A —>00 L 


-1 —61/2 b 2 0 


exp 


2bog 2 (g)i 


(3.1) 


where bo = — / 3o/2(47t) 2 , b\ = — /?i/2(47t) 4 . For large g this leads to the usual behaviour 


1 


2& 0 ln JL- 

A-MS 


+ — In (2 In 

bo V 



(3.2) 


In the 1-loop case, we set 61 = 0 in Eqs. E3J, (jSH- 

Through Eqs. dyn . 0351) . (nyrni and D, r / 2 is a function g/T and g/k^g. The 
dependence on g is formally of higher order than the computation. Numerically, of course, 
there is non-vanishing dependence, as illustrated in Fig. |21 

As usual, one may choose some “optimisation” criterion which should lead to a reduced 
/x-dependence and thus reasonable convergence. We fix g opt to be the point where the 1-loop 
coupling g | has vanishing slope (“principal of minimal sensitivity”), cf. Fig. |21 and vary then 
the scale in the range g = (0.5...2.0) x g opt around this point. Results are shown in Fig. 0 
The /x-dependence indeed decreases significantly as we go to the 2-loop level. The numerical 
2-loop value is some 20% smaller than the 1-loop value. It is comforting that the 2-loop value 
is on the side to which the “error band” of the 1 -loop result points, even though it does not 
in general lie within that band. 





























N =3 T = 2 A_ 

iNf j, i 





Figure 2: A comparison of 1-loop and 2-loop values for g 2 E /T, as a function of g/T, for a 
fixed T/A^g = 2.0 and Nf = 0, 2,3. 


N f = 0 N f = 2 N f = 3 



Figure 3: The 1-loop and 2-loop values for g 2 E /T, as a function of T /A^g (solid lines). For 
each T the scale jl has been fixed to the “principal of minimal sensitivity” point fa op t following 
from the 1-loop expression, and varied then in the range Jl = (0.5...2.0) x Jl opt (the grey bands). 


4. Spatial string tension 


The computations in the previous sections can be given a “phenomenological” application, 
by considering lattice measurements of the so-called spatial string tension. The spatial string 
tension is obtained from a rectangular Wilson loop W S (R\, R 2 ) in the {x\, X 2 )-plane, of size 
R\ x R 2 . The potential V S (R .\) is defined through 

V a (Ri) = - lim -i-In W S (R 1 ,R 2 ), (4.1) 

R 2 —KX) 1\‘2 

and the spatial string tension a s from the asymptotic behaviour of the potential, 




lim 

Ri —>00 


Vs(Ri) 

Ri 


(4.2) 
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Since a s has the dimensionality GeV 2 , it is often expressed m as the combination 



(4.3) 


where <)> is a (decreasing) dimensionless function, and T c is the critical temperature of the 
deconfinement phase transition. 

We now turn to how the result for that we have obtained in this paper, combined with 
other ingredients, allow us to obtain an independent prediction for the spatial string tension. 


4.1. Three-dimensional prediction 

The very same observable as in Eq. Oh exists also in 3d SU(3) gauge theory, or “Magne¬ 
tostatic QCD” (MQCD). Since the gauge coupling g ^ of MQCD is dimensionful, a s must 
have the form a s = c X g^, where c is a numerical proportionality constant. It has been 
determined with lattice Monte Carlo methods most recently in Ref. m where, after the 
continuum extrapolation, it was expressed as 


= 0.553(1) . 
9m 


(4.4) 


In order to compare Eqs. (EH, ((231), we need a relation between T and g^- In the previous 
section, we obtained a relation between T and g^. The relation between g ^ and g ^ is also 
known, up to 2-loop order 1251 : 2 


9m = 9e 



1 9e Ca 
48 7rmE 


17 / 9e^a \ 2 
4608 V 7T7TIE / 


(4.5) 


where the 1-loop part was determined already in Ref. 125] . 

It is worth stressing that the corrections in Eq. 631) are in practice extremely small, even 
for values of tue /corresponding to temperatures very close to the critical one. (For N c = 3 
and Nf = 0, (ute/^e) 2 ~ 0-32 log 10 (T/Ajjg)-|-0.29.) This seems by no means obvious a priori , 
given the observed slow convergence in the case of the vacuum energy density of EQCD ma¬ 
in view of this fact, however, we can safely ignore all higher loop corrections in Eq. 63d- 

Another source of errors in going from EQCD to MQCD are the higher order operators that 
have been truncated from the action of MQCD. As discussed in Ref. ma, they are expected 
to contribute at the relative order 0(g^/m^), i.e. at the same order that 3-loop corrections 
enter Eq. 63d- From this consideration, one might expect them to again be numerically 
negligible. In principle one could avoid this assumption, however: the ratio y/os /has been 
estimated in Ref. UK through direct numerical simulations in EQCD. Unfortunately the 
statistical and particularly the systematic errors appear to be non-vanishing (no continuum 

2 The 2-loop correction 5 <?m/s1 = — g^C a^{CaCf + 1)A^ + (6 Cf — Ca) Ag^]/384(7rmB) 2 was ignored in 
Ref. j25) . as it is of higher order according to 4d power counting and numerically insignificant. 
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extrapolation was carried out for this quantity), so that we prefer to follow the line starting 
from Eq. (TOl) in the following. Nevertheless it would be interesting to learn more about the 
importance of the higher order operators. 

Now, as we know g^/T as a function of T/A^g from Fig. £3 Eqs. (HU) and (TOl) allow us 
to obtain yTxJ/T as a function of the same variable. In order to compare with Eq. (14.31) . 
however, we still need to relate Aj^§- to T c . This problem has also been addressed with 4d 
lattice simulations, as we review in Sec. EU 


4.2. Critical temperature in “perturbative units” 


The determination of T c / A^g is a classic problem in lattice QCD. Two main lines have been 
followed, one going via the zero temperature string tension s/a, the other via the Sommer 
scale ro |38| . 

Values obtained for T c /y/a by various lattice collaborations are summarised in Ref. m , 
Table 7. Traditionally the values were around T c /y/(T = 0.630(5) EO], but Ref. [3?!. argues in 
favour of a slightly larger number in the continuum limit. Indeed the most precise estimate 
appears to come from Ref. EH, where T c /yTx = 0.646(3) is cited. Combining with A^g/y H = 
0.555(19) from Ref. |i2] . we are lead to 

= 1-16(4) . (4.6) 

-^MS 


The error is dominated by the one in A^g/ybr. 

A value for roT c , on the other hand, has been obtained in Ref. ESj: roT c = 0.7498(50). 
Combining with roA^g = 0.602(48) from Ref. EU (the value roA^g = 0.586(48) from a few 
lines below Eq. (4.11) in Ref. 05] is well within error bars), one obtains 


= L25 ( 10 ) ' 

2V ms 


(4.7) 


This is consistent, within statistical errors, with Eq. EU, if favouring a slightly larger central 
value. Again the error is dominated by the zero-temperature part, roAj^g in this case. In 
general it might be expected, though, that systematic uncertainties are better under control 
in the extraction of r o than of ybr, since the static potential needs to be computed only up 
to intermediate distances. 

Apart from going through \fa and ro, there is also a third possibility m ■ It is based 
on directly determining a (lattice) A-parameter from the scaling of a suitably defined renor¬ 
malised gauge coupling at the critical point, and converting at the end to the MS scheme. 
The value obtained is m 


At 


= 1.15(5) , 


(4.8) 


consistent with Eqs. (IPl) and (E3- 

To be conservative, we will consider the interval T c / A^g = 1.10...1.35 in the following, 
encompassing the central values as well as the error bars of Eqs. (E33)-(E!H|)- 


li 





1.2 


T /A srs =U0...1.35 



Figure 4: We compare fd lattice data for the spatial string tension, taken from Ref. \2tAj . with 
expressions obtained by combining 1-loop and 2-loop results for g 2 E together with Eg. 
and the non-perturbative value of the string tension of 3d SU(3) gauge theory, Eg. \f.f\ ). The 
upper edges of the bands correspond to T c /= 1.35, the lower edges to T c /= 1.10. 

4.3. Four-dimensional measurement 

The spatial string tension of 4d pure SU(3) gauge theory at temperatures above the critical 
one, as a function of T/T c , has been measured at N r = 8 in Ref. \TI (cf. Fig. 11). There are, 
of course, systematic uncertainties, both from the lack of a continuum extrapolation as well 
as from how the string tension is extracted by fitting to the large-distance behaviour of the 
static potential. Nevertheless, we expect that the results are in the right ballpark. 

Given the considerations in Secs. imipi we can thus compare the 3d and the 4d deter¬ 
minations of sJcif/T. The result is shown in Fig. 0J where is plotted. We observe 

a significant discrepancy at 1-loop level (as most recently pointed out in Ref. 1231 ), but a 
remarkable agreement once we go to 2-loop level. It is also noteworthy that the functional 
form of the 2-loop curve appears to match the behaviour of the lattice data down to low 
temperatures. 


5. Conclusions 

The main purpose of this paper has been the analytic computation of the 2-loop effective 
gauge coupling of QCD at finite temperatures, defined as a matching coefficient appearing in 
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the dimensionally reduced effective theory, EQCD. 3 The result is given in Eqs. C33M233- 
We have also determined a new contribution of order 0(g 6 T 4 * ) to the pressure of hot QCD; 
the information is contained in Eq. (HMD , and how it enters the pressure is explained in 
Ref. pm. 

The 2-loop correction we find is numerically substantial, some 20% of the 1-loop expression. 
This indicates that while perturbation theory is in principle still under control, if restricted 
to the parametrically hard modes p ~ 2nT only, it is important to push it to a sufficiently 
high order, in order to obtain precise results. 

Our expression for the effective gauge coupling has a direct “phenomenological” appli¬ 
cation, in that it allows for a parameter-free comparison of 3d MQCD and 4d full theory 
results for an observable called the spatial string tension. We find that the 2-loop correc¬ 
tion computed here improves the match between the two results quite significantly, down 
to temperatures very close to the critical one. A small discrepancy still remains but, given 
that no continuum extrapolation was taken in 4d lattice simulations, that the extraction of 
the spatial string tension may involve systematic uncertainties due to large subleading terms 
in the r-dependence of the spatial static potential 14 (?’) [IHj, and that there also has to be 
some room for residual 3-loop corrections, as well as improvements in the matching between 
EQCD and MQCD, we do not consider this discrepancy to be worrying. We do believe that 
the discrepancy can be decreased by improving systematically on the various ingredients that 
enter the comparison. 

These conclusions support a picture of thermal QCD according to which the parametrically 
“hard” scales, p ~ 2-kT, can be treated perturbatively, almost as soon as we are in the 
deconfined phase, will the parametrically “soft” scales, p ~ gT,g 2 T, require in general a 
non-perturbative analysis within one of the effective theories describing their dynamics. For 
the observable we considered here, in fact, even the colour-electric scale p ~ gT could be 
integrated out perturbatively, but it is known that this is in general not the case. We should 
like to stress that this conclusion is rather non-trivial, as there numerically is little hierarchy 
between the scales 2 irT, gT, g 2 T at the realistic temperatures that we have been considering. 
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3 Other “effective gauge couplings” can of course also be defined; for a recent review, see Ref. JTTJ . The 

difference is that in these cases all momentum scales influence the effective gauge coupling, so that perturbation 

theory cannot be reliably applied for its computation. 


13 




Appendix A. Expansions for master integrals 

Using the notation introduced in the text, the master integrals of Eq. flHEH) read, up to O(e): 


4(D = 


4(2) = /x 
4(3) = /X 


—2e 
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+ e 

1 


(4vr) 
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C(3) 


128vr 4 T 2 
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/xe 7 ® \ 

4ttT J 
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C(3) 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 
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